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ABSTRACT 

The  Naval  Command,  Control  and  Ocean  Surveillance 
Center  RDT&E  Division  (NRaD)  has  been  using  a  500 
MHz  Linear  Frequency  Modulated  (LFM)  radar  to 
collect  measurements  of  flying  aircraft.  These  data 
have  been  used  to  generate  high  resolution  Inverse 
Synthetic  Aperture  Radar  (ISAR)  images  of  the  targets 
[1].  Digital  Signal  Processing  (DSP)  hardware  has 
been  added  to  the  radar  and  algorithms  have  been 
implemented  to  perform  ISAR  processing  on  the  data 
in  real  time.  A  VME  bus  architecture  has  been 
developed  to  provide  a  scaleable,  flexible  platform  to 
test  and  develop  real-time  processing  software. 
Algorithms  have  been  developed  fi-om  a  system  model, 
and  processing  software  has  been  implemented  to 
perform  pulse  compression,  motion  compensation, 
polar  reformatting,  image  formation,  and  target  motion 
estimation. 

Keywords:  Radar,  Imaging;  ISAR;  Aircraft,  Linear 
Frequency  Modulation. 

1.  INTRODUCTION 

NRaD  has  a  long  history  of  generating  inverse 
synthetic  aperture  radar  images  of  ships  and  aircraft 
using  wide-hand,  coherent  radars  [2] [3].  ISAR 
processing  of  radar  data  has  many  advantages.  ISAR 
images  may  be  used  for  target  identification  or  to 
determine  target  configuration.  The  processing  also 
provides  coherent  integration  gain  that  will  enhance 
target  detection.  ISAR  processing  has  test  range 
applications  performing  diagnostic  measurements  of 
targets  and  determining  specification  compliance. 

Advances  in  digital  signal  processing  (DSP) 
technology  .  make  it  possible  to  perfoim  ISAR 
processing  in  a  computer  processor  that  is  embedded  in 
the  radar  system  and  generates  ISAR  images  at  the  rate 
that  the  radar  signals  are  being  received.  This  means 
that  a  radar  operator  is  able  to  view  target  images  as 
the  target  is  being  illuminated  by  the  radar. 


The  NRaD  Radar  Branch  is  developing  a  test  bed 
system  to  aid  in  the  development  of  the  architecture 
and  algorithms  needed  to  perform  real-time  ISAR 
processing.  This  system  is  described,  and  initial  results 
are  presented. 

2.  DESCRIPTION  OF  LFM  RADAR 

The  NRaD  LFM  radar  is  a  wide  band,  coherent,  high 
range  resolution  radar  that  operates  at  X-band.  It 
generates  a  linear  frequency  modulated  waveform  with 
500  MHz  bandwidth  and  a  600  ns  pulse  width.  Time- 
domain  sampling  of  the  received  signal  and  digital 
pulse  compression  yield  a  range  resolution  of  0.5 
meters  over  a  range  window  of  up  to  1200  m.  Table  1 
summarizes  the  radar  parameters. 


LFM  Radar  Parameter  Summary 

Bandwidth: 

500  MHz 

Pulse  Width: 

600  ns. 

Center  Frequency,  fc: 

9.25  GHz 

Transmit  Power,  Pt: 

1000  W  (peak) 

Pulse  Frequency,  PRF: 

50  -  2000  Hz 

Receiver  Noise  Figure,  F^: 

6dB 

Sample  Frequency,  fs: 

500  MHz 

Digitizer  Dynamic  Range: 

12  bits  ( 72  dB  ) 

Spurious  Free  Dynamic  Range: 

30  dB 

Range  Resolution: 

0.5  m.  (3  ns.) 

Table  1. 


The  radar  is  located  on  the  west  side  of  Point  Loma, 
overlooking  the  Pacific  Ocean,  at  approximately  120 
feet  above  sea  level.  This  allows  the  dynamic 
measurement  of  aircraft  at  ranges  of  1  to  5  miles  and 
altitudes  of  200  to  2800  feet.  Ship  targets  can  also  be 
measured  at  ranges  of  3  to  8  nmi.  Reference  2  contains 
a  summary  of  the  characteristics  of  the  measurement 
site. 

2.1  RADAR  DATA  COLLECTION 

The  LFM  Radar  uses  a  surface  acoustic  wave  (SAW) 
filter  to  generate  a  “chirp”  signal  that  sweeps  from 


1050  MHz  to  1550  MHz  in  600  ns.  This  is  mixed  to 
an  X-band  center  frequency  of  9.25  GHz  and 
transmitted.  The  received  signal  is  mixed  back  to  a 
center  frequency  of  1300  MHz.  The  signal  is  then 
quadrature  ph^se  detected  by  mixing  it  with  the  1300 
MHz  IF  reference  signal.  This  yields  an  in-phase  (I) 
and  a  quadrature-phase  (Q)  signal,  each  of  which  has  a 
250  MHz  bandwidth.  This  process  of  baseband 
conversion  allows  for  the  minimum  sampling  frequency 
while  still  achieving  the  Nyquist  rate  [4]. 

The  in-phase  and  quadrature-phase  detected  signals  are 
sampled  and  digitized  at  the  Nyquist  rate  of  500  MHz 
per  channel.  The  sample  clock  is  coherent  with  the 
master  oscillator,  and  the  PRF  and  range  trigger  are 
integral  numbers  of  sample  periods,  so  the  time  and 
phase  delay  from  the  transmitted  pulse  to  each  sample 
are  known.  When  triggered,  the  digitizer  can  capture 
4096  I-Q  sample  pairs.  This  corresponds  to  an  8  ms 
time  window  (1200  m.  range  window).  This  window 
size  is  selectable;  however,  it  must  exceed  the 
transmitted  pulse  length  plus  the  delay  extent  of  the 
target.  The  digitized  waveform  is  transferred  into  the 
memory  of  a  microprocessor-controlled  data  collection 
system. 

3.  ISAR  THEORY 

An  ISAR  image  is  an  estimate  of  the  radar  reflectivity 
of  a  target  with  respect  to  a  spatial  coordinate  or 
coordinates.  A  one-dimensional  image  where  the 
coordinate  lies  along  the  line-of-sight  from  the  radar  to 
the  target  is  called  a  range  profile.  Typically  in  the 
case  of  two-dimensional  images  of  aircraft,  two 
orthogonal  coordinates  are  assumed  to  lie  in  a  (nearly) 
horizontal  plane  that  contains  the  target  and  the  radar. 
The  origin  is  located  at  an  arbitrary  point  on  the 
aircraft. 

A  target  is  typically  modeled  as  a  collection  of  small, 
“point  scatterers.”  Each  scatterer  reflects  the 
transmitted  signal  scaled  in  amplitude  and  delayed  by 
2RJc,  where  R,  is  the  range  to  the  i^  scatterer.  The 
received  signal  is  the  sum  of  the  reflections  of  each 
scatterer: 

s(l)  =  I pit or 

where  p(t)  is  the  transmitted  signal,  P(co)  is  its  Fourier 
Transform,  and  a  is  the  “radar  cross-section”  of  the 
scatterer  [3]  [4]. 


3 . 1  ONE-DIMENSIONAL  IMAGE  MODEL 

A  target  function  can  be  defined  to  represent  Vcr  as  a 
function  of  some  spatial  coordinate  system.  Let  h(R-Ro) 
be  the  target  function  where  Rq  is  the  range  to  an 
arbitrary  fixed  point  on  the  target.  Then, 

5(6;)=  f  ■  ' . , 

If  Ro  is  much  larger  than  the  target  extent,  then 
S{(o)  = 

^  -DO 

= /c) 

Thus,  when  the  term  which  represents  the  effects  of 
target  range  is  removed,  the  remaining  target  response 
is  the  convolution  of  the  transmitted  signal  with  the 
target  function:  s(t)  -  p(t)  *  h(t  •  c/2).  The  one¬ 
dimensional  imaging  process  involves  inverting  this 
relation  to  estimate  h(r)  from  s(t).  This  is  a 
deconvolution  operation  [4]. 

3.2  TWO-DIMENSIONAL  IMAGE  MODEL 

The  geometry  for  the  two-dimensional  case  is  shown  in 
figure  1.  The  target  function,  h{x,y),  has  a  two- 
dimensional  spatial  Fourier  Transform  of  H{kx  If 
the  target  is  illuminated  by  a  plane  wave,  the  R  in  the 
above  equations  can  replaced  with  Rq-x  cos0  -  y  sin9. 


wave  front 

Figure  1.  Two-Dimensional  Image  Geometry 

Integrating  over  the  target  yields  a  result  similar  to 
one-dimensional  case. 

Seiw)  =  P{co)H{=2^cose,^s\ne) 

Solving  for  co  and  0  yields 
0)  =  jc^yjkl  -H  kj,  0  =  i2iX\~\ky  I  kf)  . 


This  shows  that  ©  and  0  are  polar  coordinates  in  the 
target  function’s  spatial  frequency  domain.  Again,  the 
imaging  problem  involves  inverting  the  above  equation 
to  obtain  an  estimate  of  h{x,y). 

4.  ALGORITHMS 

For  each  received  pulse,  a  vector  of  N  digitized  I  and 
Q  samples  is  collected.  The  sampled  representation  of 
the  complex  received  signal,  s(t),  is  -  jQ^,.  When 
a  Discreet  Fourier  Transform  (DFT)  is  applied,  the 
result  is  a  sampled  representation  of  S(o)).  Because  of 
the  baseband  conversion  process,  these  samples, 
correspond  to  the  radar  frequencies  from  fc  -fs  12  to 
fc‘^fsl2,  in  steps  of/^/N. 

Pulsed  radars  illuminate  the  target  (and  collect 
received  data)  periodically  over  time.  However,  this 
time  scale  is  a  much  larger  than  the  one  used  to 
characterize  the  transmitted  and  receive  signals.  It  is 
convenient  to  define  a  slow-time  variable,  u.  Functions 
of  w,  such  as  target  range  and  incidence  angle  (0)  can 
be  considered  constant  with  respect  to  fast-time,  t.  A 
frame  of  data  consisting  of  M  pulses,  each  having  N 
samples  will  form  an  MxN  matrix  where  each  row 
contains  the  successive  samples  of  a  pulse,  separated  in 
time  by  2  ns.,  and  the  rows  contain  successive  pulses, 
separated  by  the  pulse  repetition  interval  (PRJ)  of  5  to 
100  milliseconds. 

4.1  PULSE  COMPRESSION 

The  inversion  (imaging)  process  for  either  model 
above  involves  evaluating  the  equation, 

H((o)=  P(o))-’S(a)). 

This  is  most  easily  implemented  in  the  frequency 
domain.  P(o))'^  is  realizable  only  if  the  transmitted 
signal  is  infmite  in  bandwidth.  It  is  approximated  by  a 
function,  G(o)),  that  is  non-zero  over  the  band  where 
P(co)  is  non-zero.  When  G(o))  =  P(co)*  it  is  a  matched 
filter.  G(co)  can  also  be  a  pulse  compression  filter. 
When  applied  to  the  transmitted  pulse,  this  filter 
should  result  *in  an  “ideal”  compressed  pulse  with  a 
linear  phase  spectrum  and  an  amplitude  spectrum  equal 
to  a  frequency  weighting  function  chosen  to  meet 
sidelobe  requirements  [5]. 

The  pulse  compression  filter  is  derived  from  samples 
of  the  transmitted  pulses  that  are  routed  through  the 
receiver  and  collected  in  the  same  manner  as  s(t).  Pulse 
compression  is  performed  by  Fast  Fourier  Transform 
(FFT)  of  each  received  pulse,  multiplication  by  the 


compression  filter  function,  and  then  inverse  FFT.  The 
resulting  range  profile  is  the  target  function  convolved 
with  p(t)'^g(t).  Thus  the  shape  of  this  function 
determines  the  range  resolution  of  the  image.  For  a 
bandpass  signal  with  bandwidth  (3,  the  resolution  is 
approximately  c/2p.  A  trade-off  between  resolution 
and  maximum  sidelobe  level  can  be  set  by  the  choice 
of  amplitude  weighting  function. 

4.2  MOTION  COMPENSATION 

The  next  step  in  image  generation  is  motion 
compensation,  the  process  of  removing  the  effects  of 
the  target’s  radial  motion.  For  each  pulse,  the  range  to 
the  origin  point  on  the  target  contributes  the  term 

to  the  spectrum  of  the  received  signal.  In 
principle  this  can  easily  be  removed  if  the  range  in 
known  for  each  pulse.  The  phase,  however,  is  sensitive 
to  relative  range  errors  of  less  that  A,/4,  (  <  1  cm  at  X- 
band),  so  range  must  be  accurately  estimated.  The 
range,  R(u),  is  modeled  as  a  third  order  polynomial 
(i.e.,  range,  velocity,  acceleration,  and  jerk)  over  an 
image  frame.  The  compensation  is  performed  while  the 
data  is  in  the  frequency  domain. 

4.3  RANGE-DOPPLER  IMAGE 

After  motion  compensation,  any  amplitude  or  phase 
modulation  due  to  radial  motion  of  the  target  will  have 
been  removed.  However,  if  the  target  is  rotating,  the 
relative  range  of  scatterers  that  are  offset  in  cross-range 
from  the  origin  will  change  over  successive  pulses.  The 
change  in  range  will  be  2A„;sin(A0),  where  is  the 
cross-range  distance  to  the  scatterer.  If  the  angle 
change  is  small,  and  the  rotation  rate,  0',  is  constant, 
then  the  response  of  that  scatterer  will  have  a  constant 
frequency  Doppler  modulation. 

h{vi)  = 

Range-Doppler  images  are  generated  at  this  stage  by 
taking  the  Fourier  Transform  with  respect  to  u,  (i.e.  an 
FFT  by  column). 

H{kJ  =  h{Qi)5{k^  +  2o)k^,  /  c) 

This  results  in  the  response  of  the  scatterer  being 
shifted  in  the  ki  or  cross-range  direction.  The  resulting 
image  shows  target  reflectivity  as  a  function  of  range  in 
one  dimension,  and  Doppler  frequency  in  the  other 
dimension.  Scatterers  rotating  towards  the  radar  have 
positive  doppler  frequency,  and  those  rotating  away 
have  negative  frequencies. 


The  range-doppler  image  has  several  uses.  When  the 
assumptions  made  above  are  valid,  H(r,  XkJ26)  is  a 
good  approximation  of  h(x,y),  where  x  is  aligned  with 
the  radar  line-of-sight.  In  addition,  ^^does  not  have  to 
be  known  to  generate  the  image.  If  the  rotation  rate  {O') 
is  not  known,  however,  then  cross-range  dimensions 
cannot  be  calculated.  This  reveals  an  ambiguity  that  is 
inherent  in  the  two-dimensional  imaging  geometry. 

Since  processing  over  u  (slow-time)  involves 
successive  radar  pulses,  the  pulse  frequency  (PRF)  of 
the  radar  is  the  sampling  frequency,  and  thus  it 
determines  the  maximum  Doppler  frequency  that  can 
be  resolved  without  aliasing.  A  range-doppler  image 
allows  this  aliasing  to  be  directly  observed.  In  addition, 
errors  in  estimating  the  target’s  velocity  cause  an  offset 
of  the  entire  image  in  the  i,,  dimension.  This  effect  can 
be  used  as  an  aid  in  refining  motion  estimates. 

4.4  TWO  DIMENSIONAL  IMAGING 

Several  factors  contribute  to  errors  in  the  range- 
doppler  image.  As  the  change  in  angle  increases,  non¬ 
linear  components  of  scatterer  motion  causes  phase 
error.  In  addition,  if  the  change  in  range  of  a  scatterer 
exceeds  the  range  resolution,  the  scatterer’s  response  is 
smeared.  The  image  becomes  blurred  as  the  distance 
from  the  origin  increases.  Compensating  for  these 
effects  is  accomplished  by  inverting  the  two- 
dimensional  image  model. 

Pulse  compression  and  motion  compensation  of  a 
frame  of  data  results  in  a  matrix,  , consisting  of 

samples  of  H(u,co).  If  the  target’s  rotation,  0(u),  is 
known  or  can  be  estimated,  then  an  angle  can  be 
assigned  to  each  response,  resulting  in  H(6,q)).  The 
data  then  represents  the  target’s  reflectivity  as  a 
flmction  of  frequency  and  rotation  angle.  Data  in  this 
form  are  said  to  be  polar-formatted,  because  frequency 
and  rotation  angle  form  polar  coordinates  in  the  spatial 
frequency  domain,  H(kj,,ky).  To  generate  a  focused 
ISAR  image,  the  data  must  be  polar-reformatted.  The 
FFT  algorithm  requires  evenly  spaced  data  samples, 
and  polar-formatted  data  are  normally  not.  Therefore, 
the  data  must  be  resampled  at  points  on  an  evenly 
spaced  rectangular  grid  in  the  kj, ,  ky  space.  A  two- 
dimensional  inverse  FFT  is  then  applied  to  recover 
Hx,y)  [2],  '  . 

In  the  resampling  technique  currently  used  for  LFM 
data,  each  row  of  frequencies  is  interpolated  to  find 
evenly  spaced  points  in  down-range  spatial  frequency. 
Then  each  column  of  the  new  data  is  interpolated  to 


find  evenly  spaced  points  in  cross-range  spatial 
frequency  [3]. 

5.  SYSTEM  ARCHITECTURE 

The  NRaD  Real-Time  ISAR  System  is  based  around  a 
VME  bus  [6].  This  allows  various  modular  system 
components  such  as  processors,  memor>',  and  interface 
boards  to  be  added.  The  system  configuration  is  shown 
in  figure  2. 


Figure  2.  NRaD  Real-Time  ISAR  System 


Digitized  samples  of  the  received  pulse  are  transferred 
by  the  Memory  Interface  card  to  the  Global  Memory 
card  via  a  dedicated  VSB  bus.  This  can  be  done  in 
parallel  with  transfers  over  the  VME  bus.  The  time 
synchronization  board  generates  time  tags  for  each 
pulse.  The  CRT  card  is  used  to  drive  a  touch-screen 
status  and  control  panel.  The  Radar  Interface  card  is 
used  to  read  parameters  from  the  radar  and  provide 
control  outputs.  The  Real-Time  Control  Processor 
coordinates  system  activity  through  the  above  I/O 
cards,  manages  data  buffers,  and  outputs  the  data  to  8- 
mm  tape  for  later  analysis. 

The  Digital  Signal  Processor  (DSP)  card  performs  the 
ISAR  processing  algorithms.  It  operates  on  data 
samples  stored  in  global  memory,  and  puts  results  back 
into  global  memory.  The  card  contains  two  ADSP- 


21020  processors  operating  at  a  33-MHz  clock 
frequency.  This  results  in  a  peak  processing  rate  of  67 
million  floating  point  operations  per  second  on  each 
processor. 

The  Display  Processor  generates  waveform  and  image 
displays  from  the  processed  data  in  global  memory.  It 
is  a  PC  compatible  computer  that  is  also  used  for 
software  development  and  debugging.  Because  the 
system  is  based  on  a  bus  architecture,  its  capabilities 
can  be  expanded  by  adding  additional  DSP  or  display 
processor  cards  in  the  spare  slots. 

5.1  DATAFLOW 

Figure  3  shows  the  flow  of  data  through  the  Real-Time 
ISAR  processing  algorithms. 
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point  for  the  image.  The  range  delay  of  this  point  is  fed 
into  the  motion  estimation  filter. 

A  sequence  of  range  profiles  can  also  be  processed  to 
form  a  range-doppler  image.  This  is  the  least 
computationally  intense  image,  and  is  useful  for  target 
acquisition  and  tracking.  When  range-doppler  images 
are  generated,  the  doppler  frequency  of  the  prominent 
scatterer  is  determined  and  used  to  calculate  a  velocity 
correction  to  the  target  motion  estimate. 

Next  the  target  spectrums  are  collected  in  a  frame 
buffer  to  form  Em,n  ♦  A  rotation  estimate  is  used  to 
calculate  the  interpolation  points.  Currently  the 
resampling  process  does  not  attempt  to  rotate  the  x  and 
y  axes  to  an  orientation  aligned  with  line  of  flight. 
Therefore  the  image  still  shows  range  on  one  axis  and 
cross-range  on  the  other.  The  resampling  also  assumes 
a  constant  rotation  rate  as  well.  A  2-D  IFFT  of  the 
interpolated  values  results  in  an  estimate  of  h(x,y).  The 
displayed  image  is  the  magnitude  of  ,  scaled  in  dB. 

6.  RESULTS 

Currently  the  algorithms  described  above  are  being 
tested,  and  the  software  to  implement  them  is  being 
transferred  from  a  workstation  to  the  digital  signal 
processors.  The  data  being  processed  are  raw  samples 
that  were  recorded  on  8-mm  tape  by  the  LFM  radar. 
These  samples  were  collected  from  stationary 
calibration  targets,  small  boats,  commercial  airliners, 
and  military  jets.  The  use  of  pre-recorded  data  allows 
testing  without  requiring  the  radar  or  targets  to  be 
available. 

Typically  the  digitizer  collects  1024  /  and  1024  Q 
samples  per  pulse  at  500  MHz.  This  results  in  a  2  ps 
sample  time,  which  assures  that  the  entire  return  signal 
is  sampled,  even  when  only  coarse  tracking  is  being 
used.  Transferring  this  data  to  global  memory  takes 
about  1  ms.  Transferring  it  to  the  DSP  memory  also 
requires  1  ms.  These  transfers  are  done  by  dedicated 
hardware,  and  occur  in  parallel  with  the  signal 
processing. 


Figure  3.  Real-time  ISAR  Data  Flow 

I  and  Q  samples  of  the  received  signal  are  read  from 
global  memory  and  combined  into  a  complex  vector, 
The  vector  is  Fast  Fourier  Transformed  to  S„  and 
multiplied  by  the  pulse  compression  and  motion 
compensation  vectots,  resulting  in  the  target  Spectrum 
A  range  profile,  h„ ,  is  generated  by  IFFT  and  saved  in 
a  buffer  for  display  and  range  tracking  purposes.  A 
prominent  scatterer  is  tracked  and  used  as  the  origin 


The  first  DSP  performs  the  following  functions: 

1)  convert  integer  samples  to  complex,  (190  ps), 

2)  Fourier  Transform,  (590  ps), 

3)  calculate  range  phase  vector,  (125  ps), 

4)  evaluate  motion  compensation  exp.,  (2,920  ps), 

5)  multiply  compression  filter  and  motion 
compensation  vector,  (250  ps), 

6)  inverse  Fourier  Transform,  (800  ps), 

7)  square  vector,  (65  ps), 


8)  find  peak  index,  (100  fis), 

9)  and  copy  data  to  buffers,  (125  ps). 

This  requires  a  total  processing  time  of  about  5.2 
milliseconds.  A  design  goal  is  to  process  pulses  at  200 
Hz  that  would  allow  less  than  5  ms  of  processing  per 
pulse.  A  significant  portion  of  the  time  is  taken 
evaluating  the  exponential  in  the  motion  compensation 
term.  At  the  end  of  each  pulse  processing  cycle,  the 
target  spectrum  and  range  profile  data  reside  in 
memory  buffers  which  are  accessible  to  the  second 
DSP.  After  a  frame  of  data  is  accumulated,  it  performs 
the  next  steps  in  the  algorithm  while  the  first  processor 
fills  a  second  buffer. 

A  typical  frame  size  for  commercial  aircraft  is  256 
pulses,  just  over  1  second.  A  256  point  FFT  requires 
130  ps.  With  copying  and  computing  power,  a  256  x 
256  range-doppler  image  can  be  generated  in  45 
milliseconds.  If  the  target  is  small  this  can  be  reduced 
further.  The  image  data  must  be  scaled  to  a  byte  value 
and  v^itten  to  the  memory  of  a  display  card  in  order  to 
be  viewed.  This  still  allows  an  image  update  rate  of  10 
Hz. 

The  two-dimensional  inverse  FFT  required  for  the 
focused  image  requires  about  100  ms.  In  addition  the 
polar  sampling  must  be  performed.  Since  sampled 
range  is  generally  8  or  more  times  longer  than  the 
target,  the  data  in  the  co  domain  (rows)  is  over¬ 
sampled.  Thus  a  linear  interpolation  between  adjacent 
point  will  yield  good  results.  For  each  pulse,  an  angle, 
0m  is  assigned.  The  kx  sample  for  each  frequency  o)n  is 
taken  at  o)=a)n/cos0m. 

The  ky  samples  are  found  by  interpolating  along  the  0 
dimension  (columns).  The  samples  are  taken  at  a 
spacing  of  A9=  A/^cOo/cOn.  Since  the  data  is  not  over¬ 
sampled  in  the  0  domain,  this  interpolation  is 
performed  in  the  workstation  environment  using  an 
FFT,  zero  padding,  and  an  IFFT.  This  approach  is  too 
computationally  intensive  for  a  real-time  environment. 
In  addition,  evaluating  the  polar  to  rectangular 
mapping  in  cases  where  rotation  is  not  uniform  poses 
an  additional  challenge. 

*  7.  CONCLUSIONS 

There  are  definite  advantages  to  a  real-time  imaging 
capability  in  high  resolution,  coherent  radars.  The 
architecture  developed  at  NRaD  shows  promise  in 
providing  this.  Additional  work  is  required  in 
implementing  the  classical  imaging  equations  in  a 
computationally  efficient  manner.  Two  areas  that  show 


room  for  improvement  are  the  evaluation  of  the  phase 
exponential  used  for  motion  compensation  and 
interpolation  for  polar  resampling. 
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